High impact of spin-nematic order on the lattice domains 
in thin films of iron-based superconductors 
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The spin-nematic scenario proposed for the structural instability of iron-based superconductors 
is studied theoretically within a Ginzburg-Landau approach in epitaxial thin films. The transition 
temperature and the domains that appear due the interface with the substrate are computed analyt- 
ically by following a variatonal procedure. It is predicted a cross-over thickness h that separates the 
Kittel regime, in which the size of the domains goes as the square root of the film thickness I, from 
a new ultrathin regime in which the system tends to reach the single-domain state as I — > 0. The 
experimental observation of this cross-over is therefore proposed as a direct probe of spin nematiciy. 



The discovery of superconductivity in ferropnictides 
has triggered a huge interest on the properties of these 
materials [1]. It is believed that magnetism plays a 
nontrivial role not only in the superconducting proper- 
ties, but also in the structural features of these systems. 
Specifically, it has been suggested that the lattice distor- 
tion that precedes their antiferromagnetic transition is 
driven by a spin-nematic ordering of the electronic sub- 
system [2] (see Fig. 1). At the conceptual level, the 
establishment of spin nematicity is rather nontrivial. In 
the 2D Heisenberg model, for example, it is viable only 
in a certain range of parameters such that the magnetic 
ground state is determined by the minimization of fluc- 
tuations, which makes it possible the emergence of the 
spin-nematic order parameter that further can survive 
above the antiferromagnetic transition temperature [3]. 
Electron nematicity can also arise from an unequal pop- 
ulation of orbitals with different symmetries due to spin- 
orbital physics [4]. In any case, the absence of direct 
coupling with external fields makes this scenario quite 
elusive to experimental verification. Nevertheless, the 
belief in spin-nematicity has taken roots in the high-T c 
community and has motivated a flurry of studies on the 
iron-based systems. So far most of the attempts have 
targeted its manifestations in the electronic or magnetic 
properties with different degree of finality [5]. 

If the spin-nematic scenario is appropriate to describe 
the lattice distortion in ferropnictides, the nematic order 
has to manifest in other properties of the lattice as well. 
The strong response of the lattice to superconductivity 
and its softening observed in ultrasound experiments, for 
example, have recently been interpreted according to this 
picture in Refs. [6] and [7] respectively (see also [8] for 
a structural study). This interpretation, however, does 
not completely rule out the more standard scenario in 
which the lattice softens by itself. In this Letter, we indi- 
cate an experimental situation in which the spin-nematic 
mechanism for the structural distortion can be clearly 
distinguished from a conventional lattice instability. In 
particular, we show that in epitaxial thin films the spin- 
nematic mechanism is expected to have a high impact on 
the structural domains below a certain film thickness. 



In epitaxial thin films, the interface with the sub- 
strate generates stresses that have been proven effective 
to increase the superconducting transition temperature 
in FeTe [9]. These stresses are known to be responsible 
for the appearance of different structural domains, whose 
size is expected to follow the Kittel I 1 / 2 law as a func- 
tion of the film thickness I [10]. This behavior in fact has 
been confirmed experimentally in conventional ferroelas- 
tics over several orders of magnitude in film thickness 
[11]. The first exception to this rule has been recently re- 
ported for the 90° domains observed in PbTiC>3 /DyScC>3 
[12] . In this case, lattice misfit plays a key role and devia- 
tions from the Kittel Z 1 / 2 law were predicted long ago [13] . 
In our case, however, lattice misfit does not break the 
initial symmetry and therefore is expected to be unim- 
portant as discussed in [14]. Thus, within a pure elastic 
scenario as described in [15, 16], the Kittel I 1 / 2 law is ex- 
pected to be valid for the tetragonal O monoclinic insta- 
bility in the iron-based superconductors. In the following 
we show that, in contrast, if the spin-nematic ordering is 
behind this transition, this law is not obeyed in the ul- 
trathin limit. Instead, the size of the monoclinic domains 
becomes ~ I 2 /I where Z* is determined by spin-nematic 
parameters (see below). This quantity Z* also defines the 
corresponding cross-over thickness, which is experimen- 
tally accessible as we show below. 
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FIG. 1: Schematic representation of the spin-nematic sce- 
nario postulated for the iron-based superconductors. In the 
intermediate (spin-nematic) phase, there is a long-range cor- 
relation between the spin fluctuations in the Fe sublattices 
that breaks the tetragonal symmetry. The establishment of 
this spin-nematic phase from the paramagnetic state at T n is 
assumed to be behind the tetragonal o monoclinic structural 
transition in these systems. 
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The reason for the appearance of ferroelastic domains 
is the minimization of the stresses generated at the in- 
terfaces. This can be described by the Navier-Cauchy 
(or elastostatic) equations djO^ = in the absence of 
external forces, where cry is the stress tensor [17]. For 
an isotropic elastic medium this tensor is given by the 
constitutive equation <jy = XSijSkk + 2/X£y, where is 
the strain tensor and A and fi are Lame parameters. The 
strain tensor, in its turn, can be expressed in terms of 
the displacement vector u as = h(djtii + diUj), what 
will be convenient in our further analysis. Such a consti- 
tutive equation can be obtained as er^ = - ^ g £ EdV from 
the elastic free energy 

FE = \e 2 kk + ^l. (1) 

Within a Ginzburg-Landau approach, the essential in- 
gredient to describe the nematic-driven structural tran- 
sition is the free energy 

F = Fe + F ncm + Fme-i (2) 

where 

F nem = lv 2 + \v i + C f{d i r 1 ){d jV ) (3) 

describes the spontaneous appearance of the nematic or- 
der parameter r\ and 

F M E = g {Sxy + £yx)lj (4) 

the magnetoelastic coupling that eventually links the 
nematic ordering with the lattice distortion [18]. The 
coefficient a is assumed to vary with temperature as 
a = a! {T — T°) while the rest of parameters are taken 
as positive constants. Microscopically, r\ can be under- 
stood as the result of non-zero long-range correlations 
between the magnetic degrees of the Fe sublattices, say 
mi and m.2, while (mi) = (1x12) = 0. The establishment 
of spin nematicity is expected to be a 2D phenomenon 
taking place within the Fe planes with a rather weak in- 
terplane coupling. According to this, the gradient term 
in Eq. (3) can be simplified to c [(d x r]) 2 + (d y r]) 2 ] . Here 
and hereafter we assume that the xy-plane of the film 
is parallel to the Fe planes and we refer the axes to the 
crystallographic unit cell (with more than one Fe atom) 
where the structural distortion is monoclinic. 

As we have mentioned, the presence of the substrate 
makes the appearance of nematicity not necessarily uni- 
form in the film, which automatically implies the possi- 
bility of having different elastic domains since these fol- 
low the nematic ones. In fact, the nematic state has 
to be compatible with the absence of normal forces at 
free interface of the film and the vanishing of the dis- 
placements deep inside the substrate. That is, with the 
boundary conditions a xz (z = V) = o yz {z = I) = and 



u = 0. In consequence the degree of nematicity is 

z— 00 

expected to vary with the distance z to the interface with 
the substrate, and consequently in the xy plane to satisfy 
the equations of equilibrium. We assume that the film 
length and width are both much larger than any other 
characteristic length in the problem. Thus, the variation 
of nematicity in the xy plane is expected to be directly 
connected to the periodicity of the elastic domains. For 
the sake of simplicity we assume that the substrate is 
described by Eq. (1) with the same Lame parameters 
as the film, which does not change qualitatively the re- 
sults [19]. We then have elastostatic equations such that 
the nematic order parameter can be sought in the form 
1] = t/o cos(ak(z — l))e zky , which corresponds to a pattern 
of displacements u = UQCOs(ak(z — l))e lky 'k in the film 
and u = u s e kz e lky ic in the substrate as can be seen by 
direct substitution into the linearized equations of equi- 
librium [20]. This in fact gives a non-trivial solution of 
these equations when the temperature is 

T = T° + AT°( T ^-(l + a 2 )(R t ) 2 ), (5) 

where AT° = g 2 /(a'^) and h = (fie) 1 / 2 /\g\. Besides, the 
parameters k and a have to be such that 

a = cot(afcZ) (6) 

to ensure that the displacements and the normal forces 
described by this solution are both continuous at z = 0, 
which represents the matching conditions associated to 
the film/substrate interface. The actual temperature at 
which the whole system gets unstable with respect to the 
nematic phase corresponds to the highest temperature 
at which the above two conditions are satisfied, T n in 
the following. As in standard second-order transitions, 
the system is stabilized by the fourth-order term in the 
free energy (2) and the amplitude of the order parameter 
increases continuously as r]o ~ \T — T^ 1 / 2 below T n (as 
Mo and u s do as well) . This gives a nematic ground state 
that is obviously degenerate with the solution obtained 
by replacing x ■<-> y due to the symmetry of the problem. 

The transition temperature as a function of the film 
thickness is shown in Fig. 2. For large values of the film 
thickness it naturally has an asymptote that corresponds 
to the transition temperature in the free standing case 
T n .frcc = T° + AT°. In this limit, it can be said that the 
dimensions of the film are large enough to achieve a full 
relaxation of the stresses generated by the interface. The 
finite dimensions of the film in fact results in a decrease 
of the transition temperature which initially goes as T n 
(1— irl*/l)T n .f rcc . In the opposite case of small thicknesses 
the transition temperature tends to its nominal value T° . 
Here it goes as T n w T° + Q(l/l ir ) 2 in the ultrathin limit, 
where & = c/(4a'l 2 ). The crossover between these two 
regimes is determined by the quantity Z*. 
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FIG. 2: Semi-log plot of the transition temperature vs. film 
thickness (in units of I*) illustrating the expected behavior 
according to the spin-nematic scenario (black solid line). The 
blue-dashed line represents the standard behavior of conven- 
tional ferroelastics while the red one indicates a I 2 depen- 
dence. Note that the transition temperature driven by the 
spin-nematic order is restricted to the limits defined by the 
nominal transition temperature T° and the transition tem- 
perature T n .{iac of the free standing case. 



Fig. 3 displays the size of the elastic domains as a 
function of the film thickness. For thicknesses I 3> h the 
Kittel I 1 / 2 law is obeyed. In this regime the size of the 
domains is however ~ (87r) 1//2 (Z*/) 1 / 2 , so it already re- 
flects the magnetoelastic coupling through the quantity 
2*. In fact the size of the domains turns out to be similar 
to the film thickness when this is I ~ Z* , and the same 
happens with the region in which the deformation of the 
substrate is concentrated. Below this point, the system 
starts to deviate from the Kittel behavior. The reason 
is that the formation of new domains does not help to 
minimize the stresses anymore. In fact, in the ultrathin 
limit (I <C I*), these deviations are very marked and the 
domain size unusually increases as ~ nil /I with decreas- 
ing film thickness. This can be understood as due to 
the energy penalty associated with the gradients in the 
corresponding domain walls. 

These results are based on a mean-field treatment of 
the nematic transition that is expected to be valid even 
in the ultrathin limit. The reason is that the nematic ex- 
citations remain gapped at the critical point except along 
two soft lines in k space, which is due to the magnetoe- 
lastic coupling (4). This increases the effective dimen- 
sionality of the problem and makes critical fluctuations 
eventually irrelevant [16] (see also [21]). 

As we see, the quantity U defines a crossover between 
two different regimes in which the system behaves very 
differently. For I ^> 2*, on one hand, the system effec- 
tively behaves as a conventionial ferroelastic. That is, 
there is a change in the transition temperature inversely 
proportional to the film thickness and the ferroelastic 
domains that appear at this point follow the Kittel Z 1 / 2 
law. On the other hand, for I < !,, the fact that the 



FIG. 3: Log-log plot of the domain size vs. film thickness in 
units of h that illustrates the expected cross-over from the 
Kittel regime to a new quasi-single domain regime for the 
ultrathin limit. 

actual instability is driven by the nematic order parame- 
ter shows up very clearly. In this regime i) the transition 
temperature tends to its nominal value (in contrast to 
a conventional ferroelastics where it is expected to de- 
crease) and ii) the system tends to get rid the domain 
walls and the size of the domains increases by decreasing 
the film thickness. The observation of this crossover in 
iron-based superconductors would be a direct probe that 
the spin-nematic scenario is taking place. 

The precise value of 1* is quite difficult to anticipate 
due to the absence of direct measurements of the bare 
parameters that determine this quantity. Nevertheless, 
its order of magnitude can be estimated as follows. Since 
the spontaneous strains are assumed to be induced by the 
nematic order parameter, no special softness of the lattice 
itself has to be invoked, fi then can be roughly estimated 
as the atomic energy E at ~ eV per unit cell. Experi- 
mentally the spontaneous strains e = —gr\j\i saturate to 
~ 10 -3 in bulk samples. In units such that this corre- 
sponds to i] = 1, this gives the ratio \g\/n ~ 10 -3 . The 
coefficient c, in its turn, can be estimated as c ~ E n l 2 t 
taking into account that space variations of the order pa- 
rameter from 77 = 1 to —1 within the interatomic distance 
lat (~ 4A) are expected to be associated with the char- 
acteristic energy E n ~ meV of the nematic order. This 
gives 2* ~ 10 3 (-E' I i/-E l a t) 1 / 2 / a t ~ 12nm, which corresponds 
to 10 — 15 times the typical lattice parameter along the c 
direction. This conservative estimate indicates that the 
cross-over thickness 2* is accessible experimentally. 

We finally make the rather obvious remark that, by 
means of the structural domains that are expected to 
appear in epitaxial films, one can extract valuable infor- 
mation about the hypothesized spin-nematic state. At 
the mean-field level this state is described by Ginzburg- 
Landau free energy (2), which in principle can be derived 
from more microscopic approaches as the worked out in 
[22]. In Eq. (2) we have the parameter T° describing 
the nominal transition temperature that, as we have dis- 
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cussed, corresponds to the limiting value of the transition 
temperature that will be observed experimentally for the 
thinnest samples. Furthermore, the parameters a', c and 
g are associated to the quantities AT°, 1* and 9 in such 
a fashion that can all be extracted from the appropri- 
ate experimental data. The information obtained in this 
way can further be contrasted with other independent 
measurements in order to get a phenomenological char- 
acterization of the spin-nematic ordering. 

In summary, we have shown that the spin-nematic 
mechanism proposed for the structural transition in iron- 
based superconductors should reveal itself in the form 
of strong deviations from the Kittel behavior exhibited 
by standard ferroelastics. We have predicted an intrinsic 
cross-over thickness Z» below which a quasi-single domain 
regime should emerge. In this regime, the size of the 
domains is expected to be inversely proportional to the 
film thickness at the same time that the transition tem- 
perature saturates to its nominal value. We have indi- 
cated the viability of this cross-over which is connected to 
the magnetoelastic coupling in the ferropictnides. These 
findings are expected to open new routes for the study the 
spin-nematic physics, and are relevant for pseudo-proper 
ferroelastics in general. 
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